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In a lattice QCD calculation with two light dynamical Chirally Improved (CI) quarks we determine 
ground state and some excited state masses in all four A baryon channels 1/2 and 3/2 ± . We perform 
an infinite volume extrapolation and confirm the widely discussed A(1405). We also analyze the 
amount of octet-singlet mixing, which is helpful in comparing states with the quark model. 
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One of the central aims of hadron spectroscopy is to 
understand the spin-flavor-parity structure of the exci- 
tation spectra for different quantum numbers. It seems 
in particular that the nucleon and A spectra show signif- 
icant differences which, if properly understood, should 
illuminate the influence of quark mass and flavor on 
hadron structure. The lowest A(| ) mass lies below the 

Roper-like A(1600, \ + ) and the negative parity nucleon 
TV* (1535); unlike the nucleon sector it does have standard 
level ordering lying between the positive parity ground 
state and the first positive parity excitation. A baryons 
can be flavor singlets or octets or, due to the difference in 
light and strange quark masses, mixtures of both. Vari- 
ous continuum model studies discuss that mixing. This is 
the first lattice QCD analysis of the A baryons for dynam- 
ical quarks, that includes ail four, namely the J p = i - 

and | -channels. We obtain ground states compatible 
with experiment in three of those. We also study the 
infinite volume limit and give for the first time lattice 
results on singlet-octet composition for all four sectors, 
obtaining the mass of the A(1405) and confirming its fla- 
vor singlet nature. 

Lattice studies in the quenched case generally had 
problems to find the low lying A (1405). Even a study 
with two dynamical light quarks [l], Q found too large 
mass values. Only recent results Q obtained with PACS- 
CS (2 + l)-flavor dynamical quarks configurations Q 
show a level ordering which is compatible with experi- 
ment. 

We study the baryons for a set of seven ensembles with 
two dynamical Chiral Improved (CI) quarks [E|, (gj . The 
CI fcrmion action consists of several hundred terms and 
obeys the Ginsparg- Wilson relation in a truncated ap- 
proximation. The pion mass ranges from 255 to 596 
MeV, the lattice spacing lies between 0.1324 and 0.1398 
fin. The bulk of our results were obtained for lattices of 
size 16 3 x 32. For two ensembles with light pion masses 
(255 and 330 MeV) we also used 24 x 48 lattices. Thus 
m^L which controls the finite size effects is 4.08 and 5.61 
in these two situations. Further details on the action and 



our simulation, as well as results for mesons, can be found 
in 0, [|| . The strange quark is introduced as a valence 
quark and its mass fixed by the J7-mass. 

The Dirac and flavor structure of the interpolating 
fields used in our study is motivated by the quark model 
[j| [l(| , see also 11 1 . Within the relativistic quark model 
there have been many determinations of the hadron spec- 
trum, based on confining potentials and different assump- 
tions on the hyperfine interaction (see, e.g;, 12h14| ) . The 



singlet, octet and decuplet attribution [llj of the states 
has been evaluated based on such model calculations, 
e.g., in [l5| (see also the summary in [HI). 

For the A baryons we use sets of up to 24 interpolating 
fields in each quantum channel. The singlet and octet 
combinations of Table U are combined with three possi- 
ble choices of Dirac matrices (Y\, T 2 ) = (1, C75), (75, C) 
and (il,C7t75) (denoted by xii Xi and X3 f° r short) 
for J = 4 and eight combinations of Gaussian smeared 
quarks (l7l.[l8| with two smearing widths (n, w). The op- 
erator numbering is given in Table [II] All interpolators 
are projected to definite parity and all Rarita-Schwingcr 
fields (spin | interpolators in Table [TTJ) are projected to 
definite spin | using the continuum formulation of the 



J. C denotes the charge 
the time and the spatial 



Rarita-Schwinger projector [l! 
conjugation matrix, 74 and 74 
direction Dirac matrices. 

For point like quark fields, Fierz identities reduce the 
actual number of independent operators (see, e.g., [20|b 
In particular, there arc no non-vanishing point-like in- 
terpolators for A(|) and singlet A(|). We use different 
quark smearing widths such that the Fierz identities do 
not apply. Hence xi> X2 and X3 are all independent, for 
J = I all interpolators are non- vanishing. 

From the cross-correlation matrix Cik(t) = 
(Oi(t)Ok(0y) we obtain the energy levels with help of 
the variational method [U [22| ■ One solves the gener- 
alized eigenvalue problem C{t)u n {t) = X n (t)C(to)u n (t) 
in order to approximately recover the energy eigenstates 
\n). The eigenvalues A„(i) ~ exp(— E n t) allow us to 
get the energy values and the eigenvectors serve as 
fingerprints of the states, indicating their content in 
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TABLE I. Baryon interpolators: The possible choices for the 
Dirac matrices T\ 2 in the spin ^channels are discussed in the 
text. Summation convention applies; for spin | observables, 
the open Lorentz index (after spin projection) is summed after 
taking the expectation value of correlation functions. 



terms of the lattice interpolators. The quality of the 
results depends on the statistics and the provided set of 
lattice operators. The dependence on to is used to study 
the systematic error; in the final analysis we use to = 1 
(with the origin at 0). The statistical error is determined 
with single-elimination jack-knife. For the fits we use 
single exponential behavior but check the stability with 
double exponential fits; we take the correlation matrix 
for the correlated fits from the complete sample 0. 

For the extrapolation to the physical pion mass we fit 
to the leading order chiral behavior, which is linear in 
■m\. Two ensembles (at pion masses 255 MeV and 588 
MeV) suffer from a slight mistuning of the strange quark 
mass, which are therefore discarded in the extrapolation 
to the physical pion mass, whenever the effects are found 
significant. This is the case for the lowest energy levels in 
each channel (three lowest ones in h ). The quoted sys- 
tematic errors for these levels include the corresponding 
deviation and the dependence of the energy levels on the 
choice of interpolators and fit ranges for the eigenvalues. 

In the present study we are restricted to 3-quark oper- 
ators for the baryon. Note that ideally one should take 
into account also meson-baryon interpolators (see, e.g., 
the discussion in f23|). This leads to many more con- 
tributing graphs and necessitates also the inclusion of 
backtracking quark loops. The related computational 
and algorithmic effort prevented such lattice calculations 
so far, although such studies arc in progress [24j |, Due 
to sea quarks, in principle, 3-quark operators have over- 
lap with meson-baryon states as well. The corresponding 
coupling was however found to be weak in actual simu- 
lations ]3, HH . We will argue below that in particular in 
the s-wave channels we find hints of such coupling even 
for our interpolators. 



TABLE II. Baryon interpolators: Quark smearing types (n/w 
for narrow/wide) and naming convention for the interpolators 
in the different channels. The three columns for the J = \ 
interpolators refer to Xi~X3- 



(1,2,11,20,25,26,33,34,43). After extrapolation to the 
physical point our lowest energy level agrees well with 
the experimental A(1116). The systematic error esti- 
mated from different combinations of interpolators and 
fit ranges is indicated in the summary Fig. [B] Analyzing 
the eigenvectors, we find that the ground state is domi- 
nated by octet interpolators of Dirac structure \i an d X3 
in agreement with a relativistic quark model calculation 
15(). Our first excitation is dominated by singlet in- 
terpolators (first Dirac structure) matching the A(1810) 
(singlet in the quark model). The Roper-like A(1600) 
(octet in the quark model) seems to be missing. This re- 
sembles the situation in the N(^ + ) channel Q. The sec- 
ond and third excitations are again dominated by octet 
interpolators. 

Infinite volume extrapolation: We performed a vol- 
ume analysis for several sets of interpolators and various 
fit ranges. The results in finite volume and the infinite 
volume extrapolations for the ground state for specific 
interpolators are shown in Fig. [2j The extrapolation fol- 
lows the method of [26|. A stable choice is the set of 
interpolators A and t m i n = 5a. The corresponding sys- 
tematic error estimated from different interpolators and 
fit ranges is indicated in the summary Fig. H3 Our final 
result is 1126(17) (11) MeV (statistical and systematic er- 
ror), which agrees nicely with the experimental A(1116) 
mass. 

J p = ' 



Finite volume: 



J 1 



i + 



Finite volume: In Fig. [T] we show re- 
sults for the four lowest eigenstates for interpolator set 



We use different sets of in- 
terpolators and fit ranges. We stress that our basis is 
large compared to that of other studies with three types 
of Dirac structures for both singlet and octet interpola- 
tors. We can extract the lowest four energy levels, shown 
in Fig.H using interpolators (2,3,10,18,26,27,34,42). We 
find that the ground state energy level agrees well with 
the experimental A(1405). The dependence of the lev- 
els on the tuning of the strange quark mass appears to 
be sizeable, albeit it an accident of our simulation. This 
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FIG. 1. Energy levels for A spin | + in a finite box of linear 
size L m 2.2 fm. Stars denote the experimental values [l(|, 
other symbols denote results from the simulation. The full 
lines show the extrapolations linear in the m\, the broken 
curves indicate the statistical error bands. 
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FIG. 2. Systematic error of the A ^ + ground state energy lev- 
els. The levels are shown for different choices of interpolators 
and fit ranges, labelled on the horizontal axis. E.g., "A4" de- 
notes the set of interpolator "A" and a fit range for the eigen- 
values from t — 4a to the onset of noise. "A" denotes inter- 
polators (2,3,10,18,26,27,34,42), "B" denotes (3,11,18,27,34). 
For each set of interpolators and fit range, results for small to 
large lattices (spatial size 16, 24 for ensemble A66, and 12,16, 
24 for C77, for notation see are shown from left to right. 
The corresponding infinite volume limits are the rightmost 
points. 



is one reason of the large systematic error depicted in 

Fig.m 

The chosen set of interpolators is particularly suitable 
for an analysis of the content of the states. The spatial 
support of the quark fields is equivalent in all interpola- 
tors and hence does not require additional rcnormaliza- 
tion when comparing the contribution of different inter- 
polators to the eigenstate. In addition, we use several 
interpolators for each given combination of flavor and 
Dirac structures, which allows us to identify a possibly 
higher number of states with similar structure. Within 
the basis used, the ground state is dominated by sin- 



glet interpolators of all three Dirac structures. There is, 
however, a considerable mixing with octet interpolators 
(second Dirac structure) of 15-20% in ensemble A66, i.e., 
at « 255 MeV (see Fig. |4|). This mixing is expected 
to increase towards the physical point, which may com- 
plicate the functional dependence of the energy levels on 
the pion mass. The first and second excitation are both 
dominated by octet interpolators, by the second and first 
Dirac structure, respectively. 

The first and second excited energy level are both a bit 
low but compatible with the experimental A(1670) and 
A(1800). In general in the J p = \ channel one may ex- 
pect coupling to 7r£ and KN in s-wave. In 27j,|28| the ex- 
pected energy levels in finite volumes are discussed based 
on a continuum hadron exchange model. There (with 
physical hadron masses), the low- lying scattering state 
levels in the \ channel are well separated for m^L < 3. 
For the pion masses of our study, the non-interacting 
meson-baryon thresholds lie close but (except for one 
point) above the lowest energy level observed. E.g., for 
the lowest pion mass, the values are + w 1.52 
GeV, tun + tuk ~ 1-62 GeV, above the lowest level. 
This resembles the situation in the ) channel. Ear- 

lier work argued that the coupling of single baryons to 
baryon-meson channels may be too weak to lead to ob- 
servable effects jzLH^. However, in our case, in s-wave 
scattering, we cannot exclude that one or even two of the 
observed three lowest energy levels correspond to scatter- 
ing states. Note that the measured ground state energy 
level is always (except for one point) below s-wave thresh- 
old, thus supporting the A(1405) identification. 

It has been conjectured from Chiral Unitary Theory 
that the lowest state may have a double-pole [1^, [2Sfl 
and a identification strategy for lattice simulations is sug- 
gested in [30j ]. This would require asymmetric boxes or 
moving frames. 

Infinite volume extrapolation: We study the volume 
dependence of the three lowest states for different sets of 
interpolators and various fit ranges. We emphasize that 
the stability of the signals of the excitations is compa- 
rable to the ones of the ground state. The volume de- 
pendence of all three low states appears fairly flat in our 
simulation, in a few cases showing even negative finite 
volume corrections. These features are compatible with 
significant contributions of an attractive s-wave scatter- 
ing state. For interpolators (2,3,10,18,26,27,34,42) and 
imin = 5a, after infinite volume extrapolation, we show 
the result of the final extrapolation of the ground state 
energy level to the physical pion mass in Fig. [SJ The fi- 
nal result for the ground state agrees very well with the 
experimental A(1405). Both the first and the second ex- 
citation appear to be a bit low but are compatible with 
the experimental A(1670) and A(1800) (see Figs. |3] and 
[5]) and might also possibly be s-wave 7r£ and KN signals. 
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FIG. 3. Energy levels for A spin \ in a finite box of linear 
size L ~ 2.2 fm, for legend see caption of Fig. [T] Fits are 
omitted for clarity. 
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FIG. 4. Eigenvectors for A spin 1/2" ground state and first 
excitation for ensemble A66 (m n w 255 MeV). The ground 
state is dominated by singlet interpolators of all three Dirac 
structures. The first (and also the second excitation, not 
shown) is dominated by octet interpolators. Note that a con- 
siderable mixing of singlet and octet is observed (15-20% for 
the ground state). 



J p = | , Finite volume: In spin | channels, for sym- 
metric quark fields, singlet interpolators vanish exactly 
due to Fierz identities. Wc use different quark smear- 
ing widths and thus circumvent the Fierz identities con- 
structing singlet interpolators nevertheless. We derive re- 
sults for the lowest three energy levels of the variational 
analysis of interpolators (2,9,10,16). Only the ground 
state can be clearly identified and its extrapolation agrees 
with the experimental A(1890). Within the finite basis 
used, this state is dominated by octet interpolators. 

The first excitation is dominated by singlet interpola- 
tors with non-negligible octet contributions at our light- 
est pion mass (see Fig.[5])- We want to emphasize the im- 
portance of singlet interpolators for the low lying states 
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FIG. 5. Eigenvectors for A spin 3/2 + ground state and 
first excitations for ensemble A66 (m^ w 255 MeV). We em- 
phasize the domination of singlet interpolators for the first 
excitation. Such interpolators are non- vanishing only for bro- 
ken Fierz identities, which is realized by the use of different 
quark smearing widths. 



in this channel, even though those interpolators are van- 
ishing exactly for symmetric point-like quark fields. 

Infinite volume extrapolation: Within errors we do 
not observe a clear volume dependence. The final re- 
sult agrees with the experimental A(1890) mass, but with 
large uncertainty. 

J p = I , Finite volume: We choose interpolators 
(2,7,9,10,15) and find a clear gap between ground state 
and excitations. The extrapolation of the ground state 
energy level lies clearly above the experimental ground 
state A(1520) and is compatible with the A(1690). The 
excitations extrapolate close to the A(2325). 

From the eigenvectors we find that the two lowest 
states are dominated by octet, the second excitation by 
singlet interpolators. The quark model shows for A(1520) 
singlet dominance (like for its companion A(1405)). We 
conclude that we might miss a signal for the ground state 
altogether, or, alternatively, there is a strong deviation 
from the leading chiral behavior towards the physical 
point. The latter case is intriguing as it might be re- 
lated to strong coupling to an s-wave 7rE(1385) state, 
which is discussed, e.g., in [3ll. |32|. 



Infinite volume extrapolation: The volume depen- 
dence turns out to be fairly flat, in a few cases even com- 
patible with negative finite volume corrections. The final 
result in the infinite volume limit at the physical point 
again misses the experimental A(1520) and agrees with 
the A(1690) mass. 

Summary: We present a comprehensive study of spin 
i and I A baryon ground states and excitations, utilizing 
a large basis of interpolators in the variational analysis 
including differently smeared quark sources (which allows 
to sidestep the Fierz identities), three Dirac structures, 
and singlet and octet forms. Fig. [5] shows the results for 
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FIG. 6. Energy levels extrapolated to the physical pion 
mass in finite volume L « 2.2 fm (lhs) and low lying levels 
after infinite volume extrapolation (rhs). The horizontal lines 
or boxes represent experimentally known states [l6l |. where 
the box size indicates the statistical uncertainty of la. The 
statistical uncertainty of our results is indicated by error bars 
of la, second bars closely to the right indicate systematic 
errors estimated from the use of different sets of interpolators, 
fit ranges of the eigenvalues and the strange quark mass. Grey 
symbols denote a poor xVd.o.f. > 3 of the chiral fits. 



ground states and excitations for lattices of linear size 
L m 2.2 fm (lhs) and the results for the ground states in 
the infinite volume limit (rhs). Systematic errors from 
the choice of interpolators, the fit ranges of the eigenval- 
ues and the tuning of the strange quark mass have been 
investigated. In both ^ channels and in the ^ chan- 
nel we find ground states extrapolating to the experi- 
mental values, in particular we reproduce A(1405) and 
also find two low-lying excitations. In our simulation, 
A(1405) is dominated by singlet contributions, but at 
ttt-tt w 255 McV octet interpolators contribute roughly 15- 
20%, which may increase towards physical pion masses. 
The observation of A(1405) with the employed basis sug- 
gests that this state has a non-negligible singlet 3-quark 
content. The Roper-like (octet) state A(1600), on the 
other hand, may not couple to our 3-quark interpolators. 

We analyze the volume dependence and find that only 
i + 

the ground state of spin ^ shows a clear exponential 
dependence as expected for bound states. For all other 
discussed states, the volume dependence is either fairly 
flat or obscured by the statistical error. For the i + , i 
and ^ channels the remaining small deviations at the 
physical point can be easily caused by the continuum 
limit and/or dynamical strange quarks. The discrepancy 
for the | ground state might point to some effect which 
is so far not properly accounted for. 
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